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Resumen 


La frecuencia y magnitud de los eventos climáticos extremos de 
precipitación han aumentado de forma significativa en varios países del 
mundo, incluido Perú. Estos hechos causan pérdidas económicas y 
humanas, especialmente en países en vías de desarrollo. La información 
y metodologías que permitan prevenir o diseñar estrategias para 
afrontarlas son escasas o inexistentes. El objetivo de esta investigación 
fue analizar la capacidad del producto satelital IMERG (Integrated Multi- 
satellitE Retrievals) del satélite GPM (Global Precipitation Measurement) 
y datos observados a partir de estaciones meteorológicas mediante un 
enfoque mixto para estimar la distribución de lluvias extremas en la 
región del Cusco, ubicado al sur del Perú. Este enfoque mixto aprovechó 
las ventajas que ofrecen ambas fuentes de información, como es la solidez 
de los datos observados a lo largo de varios años y la resolución temporal 
horaria del producto satelital. La metodología se basó en una curva de 
crecimiento de cada región homogénea, factor de corrección y parámetros 
que estiman la función intensidad y duración para toda la región Cusco. 
Los resultados se evaluaron mediante validación cruzada entre los valores 
de precipitación diaria obtenidos del producto IMERG, enfoque mixto y 
precipitación observada para periodos de retorno de 2, 5, 10, 20, 50, 100, 


200, 500 y 1 000 años. Los resultados sugieren que la combinación de 
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datos observados de lluvia y del satélite IMERG puede ser una alternativa 


para estimar lluvias extremas en la región Cusco. 


Palabras clave: precipitaciones máximas, regiones homogéneas, GPM- 
IMERG VO6, Índice de avenida, OMM, curvas IDF. 


Abstract 


The frequency and magnitude of extreme climatic precipitation events 
have increased significantly in several countries of the world, including 
Peru. These events cause economic and human losses, especially in 
developing countries. Information and methodologies to prevent or design 
strategies to deal with them are scarce or non-existent. The aim of this 
research was to analyze the capacity of the IMERG (Integrated Multi- 
satellitE Retrievals) satellite product of the GPM (Global Precipitation 
Measurement) and observed data from meteorological stations using a 
mixed approach to estimate the distribution of extreme rainfall in Cusco 
region located in southern Perú. This mixed approach took advantage of 
both sources of information, such as the strength of the data observed 
over many years and hourly temporal resolution of the satellite product. 
The methodology was based on a growth curve for each homogeneous 
region, correction factor and parameters that estimate the intensity and 
duration function for the entire Cusco region. The results were evaluated 
by cross-validation between daily precipitation values obtained from the 
IMERG product, mixed approach and observed precipitation for return 
periods of 10, 20, 100, 500 and 1 000 years. The results suggest that the 
combination of observed and rainfall data from the IMERG satellite may 


be an alternative to estimate extreme rainfall in the Cusco region. 
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Introducción 


Las proyecciones climáticas indican que los eventos de precipitaciones 
extremas serán más intensos y frecuentes en varios lugares del mundo, 
incluyendo Sudamérica y Perú (IPCC, 2015; UNESCO, 2020; MINAM, 
2016). En un informe conjunto entre la Oficina de las Naciones Unidas 
para la Reducción del Riesgo de Desastres y el Centro de Investigación 
sobre la Epidemiología de los Desastres (CRED, 2018) se manifiesta que 
se registraron a nivel global 3 148 fenómenos naturales que provocaron 
desastres por inundaciones entre los años 1998 y 2017. Las inundaciones 
fueron el desastre de origen natural más común, representando el 43 % 
de todos los registros. Las inundaciones también están aumentando en 
frecuencia e intensidad, y se espera que continúen a la alza debido al 
cambio climático (UNESCO, 2015). 


Las precipitaciones extremas en el Perú causaron pérdidas humanas 
y económicas en los últimos años debido a las inundaciones repentinas y 
la presencia de flujos de detritos (“"huaycos”, como se denominan en Perú) 


en todo el país, siendo Cusco la región con mayor presencia de eventos 
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de huaycos, cuyo detonante suelen ser las precipitaciones intensas 
(SENAMHI, 2019). Algunos de estos impactos están relacionados con la 
fase positiva de El Niño Oscilación del Sur (ENOS) en el océano Pacífico 
(Lavado € Espinoza, 2014). Sin embargo, las precipitaciones extremas en 
el Perú están también relacionadas con el rol que ejerce el océano 
Atlántico (Espinoza et a/., 2009). Por ejemplo, se describen inundaciones 
extremas presentadas en enero-febrero de 2010 desencadenadas por 
precipitaciones extremas en la cuenca del río Vilcanota, en la región 
Cusco, al sur de los Andes peruanos (Lavado, Taylor, Sven, € Guyot, 
2012). 


Entre 2003 y 2009, las emergencias registradas en el territorio 
peruano ante la ocurrencia de fenómenos de remoción en masa como 
deslizamientos y derrumbes a causa de eventos hidrometeorológicos 
extremos alcanzaron 1 391 emergencias y 187 000 personas afectadas, 
16 000 personas damnificadas, 4 000 viviendas afectadas y más de 3 000 
viviendas destruidas (INDECI, 2012). Se estima que el 40 % del territorio 
de las regiones Cajamarca, Cusco y Huancavelica, ubicadas en los Andes 
peruanos, tendrían muy alta susceptibilidad a movimientos de remoción 
en masa a raíz de eventos hidrometeorológicos extremos (MINAM, 2015; 
Villacorta € Fidel, 2012; INGEMMET, 2010). 


Para prevenir eventos de inundación, una herramienta útil es la 
estimación de caudales de los ríos mediante modelación hidrológica de 
cuencas; no obstante, la escasez de datos de lluvia representa un 
problema. Además, si los datos requeridos son a nivel horario, la 
complejidad de obtenerlos aumenta, pues se requiere de la capacidad de 


reconocer la distribución espacial y temporal de datos de precipitación 
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(Espinosa-López, Ibáñez-Castillo, Arteaga-Ramírez, 8 Galeana-Pizaña, 
2020). En general, los datos obtenidos a partir de un pluviógrafo se 
consideran ideales para el análisis de precipitaciones extremas, pero 
debido a su cobertura espacial y tiempo de registro limitado, la 
regionalización puede resultar errónea, en especial en los Andes peruanos 
(Llauca et al., 2021) debido al papel que ejercen los Andes sobre la 
variabilidad espacial y temporal de las precipitaciones (Espinoza et al., 
2009). El enfoque tradicional del análisis de frecuencia de eventos 
extremos se basa exclusivamente en el uso de información local y tiene 
como inconveniente las escasez de series de registros históricos, así como 
la incertidumbre asociada con la calidad (Fernández 8: Lavado, 2017). De 
acuerdo con Hosking y Wallis (1997), la escasez de datos en el tiempo se 


puede atenuar con su abundancia en el espacio. 


La simulación de caudales diarios a partir de modelos de lluvia- 
escurrimiento es a menudo insatisfactoria debido a que la variabilidad 
espacial de la lluvia está pobremente representada en regiones donde los 
datos son escasos y la topografía compleja; tal es el caso de la región 
andino amazónica (Collischonn, Collischonn, 8 Tucci, 2008; Zubieta, 
Getirana, Espinoza, Lavado, € Aragon, 2017). Una fuente alternativa de 
datos diarios ante estos inconvenientes es el uso de lluvia estimada por 
sensores remotos (SENAMHI, 2016). De hecho, el empleo de datos de 
precipitaciones subdiarios puede ser potencialmente útil para simular 
caudales en cuencas andinas, donde los episodios cortos de lluvia 
convectiva son más relevantes para la variabilidad hidrológica (Zubieta et 
al., 2019). 
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En este contexto y con fines de prevención, la estimación de eventos 
de precipitación extrema máxima en 24 horas (magnitud y frecuencia) 
resulta de importancia. Por ejemplo, es útil para el diseño de estructuras 
hidráulicas, zonificación de inundaciones, y estimación económica de 
proyectos de protección contra inundaciones (Ngongondo, Xu, Tallaksen, 
Alemaw, 8: Chirwa, 2011; Noto 8: La Loggia, 2009). Por lo tanto, el interés 
suele centrarse en eventos extremos con periodos de retorno (T) mayores 
a 25 años debido a sus riesgos asociados. Para estimar T, el enfoque 
probabilístico mediante el ajuste a una distribución de función de densidad 


de probabilidad (FDP) es el más utilizado. 


De manera reciente se documenta con mayor frecuencia el uso del 
análisis regional de frecuencia basado en los L-momentos (ARF-LM) para 
precipitaciones extremas (Cassalho et al., 2017; Malekinezhad 8 Zare, 
2014; Ngongondo et a/., 2011; Yang et a/., 2010). El ARF-LM intercambia 
espacio por tiempo; proporciona un marco para la caracterización de la 
distribución de frecuencia de la precipitación máxima, y puede 
proporcionar estimaciones más precisas de cuantiles de lluvia. Por lo 
tanto, se usa el índice de avenida como una cantidad representativa de 


análisis de frecuencia regional (Hosking 8: Wallis, 1997). 


El análisis regional de frecuencia se usa en lugares donde existe 
poca o ninguna información; se propusieron varias técnicas de estimación 
regional y se aplicaron en distintas partes del mundo; sin embargo, existe 
consenso en que el análisis regional de frecuencia es confiable (Smith, 
1989; Bradley, 1998; Viglione, Laio, 8: Claps, 2007; Malekinezhad € Zare, 
2014). 
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En cuanto al uso de información de precipitación provenientes de 

productos satelitales, las primeras experiencias en investigación indican 
un mejor rendimiento del producto IMERG del satélite GPM en 
comparación con otros productos (Llauca et a/., 2021; Salas, Rosales, 
Alvarez, 8 Asención, 2020; Lavado-Casimiro et a/., 2020). Sin embargo, 
los rendimientos varían sustancialmente entre diferentes condiciones 
topográficas y climáticas, en particular en regiones dominadas por 
cambios en la orografía (Asurza, Ramos, €: Lavado, 2018). Además, aún 
existe incertidumbre acerca de la utilidad del satélite GPM-IMERG para la 
adecuada estimación de datos de lluvia sobre zonas andinas en regiones 
tropicales debido al rol que ejerce la topografía en la distribución espacial 
y temporal de la lluvia en los Andes (Manz et al., 2017; Zubieta et al., 
2017). 


En 2019, GPM actualizó sus algoritmos de datos para calibrar e 
incorporar los datos de la misión de medición de lluvias tropicales (TRMM) 
en su registro de datos (Huffman et a/., 2020). TRMM fue un satélite en 
funcionamiento de 1997 a 2015, cuyos datos se evaluaron mediante 
investigación (Mahavik 8 Tantanee, 2021). El algoritmo IMERG fusiona 
datos de las misiones TRMM y GPM, lo que brinda a los meteorólogos e 
investigadores acceso a un registro de precipitación de 20 años desde el 
1% de junio de 2000 hasta la actualidad (Tan, Huffman, Bolvin, 8 Nelkin, 
2019). 


Evaluaciones de la calidad de la información de precipitación 
provenientes de productos satelitales realizadas en la zona central de EUA 
muestran que los datos del producto IMERG Final Run están cerca de 


representar los valores observados en tierra, mientras que las 
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estimaciones de los productos IMERG Early Run y el WRF (Weather 
Research and Forecasting) presentan valores sesgados (Zhang, Lin, 8 
Bras, 2018). En China, en términos de precipitación diaria, varias 
investigaciones demostraron que el producto IMERG VO5 por general 
representa mejor las estimaciones de precipitación diaria con respecto a 
los productos TMPA-3B42V7 (Multi-Satellite Precipitation Analysis) y 
TRMM (Yuan et a/., 2019). Wang, Tang, Han, Guo y Hong (2018) llegaron 
a conclusiones similares y obtuvieron un 20 % de mejora en relación con 
los mismos productos en la misma región de análisis. Tang et al. (2020) 
y Peng et al. (2020) encontraron que el producto IMERG VO6 superó tanto 
estadística como hidrológicamente a productos satelitales como CMORPH 
(Climate Prediction Center Morphing Method); TMPA-3B42V7; CDR 
(PERSIANN-Climate Data Record), y al producto CHIRPS (Climate Hazards 
group InfraRed Precipitation with Stations). Otros estudios, como los de 
Khodadoust, Saghafian y Moazami (2017), muestran resultados similares 
al comparar el producto IMERG con el producto TRMM en Irán, incluso 
recomiendan el producto IMERG para su aplicación en países donde se 


carece de una adecuada instrumentación pluviométrica. 


En la modelación hidrológica de la cuenca del Amazonas en Perú y 
Ecuador se utilizó información de precipitación del producto GPM-IMERG 
VO3D Final Run; se encontró que para la cuenca del Ucayali —donde está 
la región Cusco— y la cuenca del Vilcanota (motivo de la aplicación de 
este estudio) tanto el producto IMERG como el producto TMPA V7 y TMPA 
RT representan de modo adecuado los caudales observados (Zubieta et 
al., 2017). 
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Para utilizar los datos de productos satelitales es necesario 
evaluarlos de forma previo en cuanto a datos de pluviómetros (Moazami 
8 Najafi, 2021). Así, se pueden usar como entrada a modelos hidrológicos 
para finalmente evaluar su desempeño en la representación de caudales. 
Asurza et al. (2018), al usar el producto IMERG para la modelación del río 
Huancané en Puno-Perú, encontraron valores del coeficiente de eficiencia 
Nash-Sutcliffe de 0.90 y sesgo de -30.1 %. Asimismo, en la cuenca del 
río Vilcanota en Cusco-Perú, Llauca et al. (2021) utilizó el producto 
satelital IMERG Early Run (debido a que IMERG demostró tener la mejor 
capacidad para representar el ciclo de precipitación diaria por encima de 
15 mm/día) para la modelación hidrológica con GR4H (Génie Rural a 4 
paramétres Horaires), con resultados alentadores en relación con 
productos como GSMaP-NRT (Global Satellite Mapping of Precipitation 
Near Real Time), CMORPH y HE (HydroEstimator), dado que presentó el 
menor sesgo en precipitación horaria en el periodo de lluvias (noviembre 


a marzo). 


Esta investigación tiene como objetivo proponer una metodología 
con enfoque mixto, el cual aprovecha la información de las estaciones 
meteorológicas y datos obtenidos por satélite para la caracterización de 
lluvias extremas en la región Cusco-Perú. La metodología propone 
complementar las series de precipitación observadas con datos obtenidos 
a partir del satélite, la cual presenta registros de periodos de corta 
duración, pero con buena cobertura espacial. De esta forma, se 
construyen curvas intensidad-duración-frecuencia (IDF) para toda la 
región Cusco, incluso en lugares donde no existen estaciones 


meteorológicas. 
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Área de estudio 


La región de Cusco está ubicada en la parte sur oriental del Perú, en las 
coordenadas -13% 31' 21.5" y -71% 58.04'. Presenta una extensión 
aproximada de 71 986 km?, con elevaciones que van de 254 a 6 268 


msam. 


Limita al norte con las selvas de Junín y Ucayali; al oeste, con la 
selva de Ayacucho y la sierra de Apurímac; al sur, con las zonas altas de 
Arequipa y Puno; y al este, con el gran llano amazónico de Madre de Dios 
(Figura 1). La temperatura máxima en la región Cusco varía de 8 a 12 *C 
sobre las zonas más altas de la cordillera, típicos del clima frío semiseco, 
principalmente en la parte centro y sur de la región Cusco, hasta máximos 
que varían de 24 a 32 “C sobre la región de selva alta y baja ubicada en 
la parte norte de la región Cusco, la cual se caracteriza por un clima cálido 
húmedo a muy húmedo. La precipitación media anual varía de 200 mm 
en zonas de alta montaña a 5 000 mm en zonas amazónicas (PACC € 
SENAMHI, 2012). 
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Figura 1. Región Cusco y estaciones meteorológicas utilizadas en el 


análisis regional de frecuencias. 


El departamento de Cusco presenta las regiones de sierra y selva, 
que en su transición varían a valles interandinos y montañas de modo que 
a mayor altitud se caracteriza por presentar lluvia y frío con humedad, así 
como clima glaciar, con hielo perenne y temperaturas muy bajas (sur del 


Cusco); a menor altitud (noreste del Cusco) se caracteriza por zonas con 
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caracteres climáticos que van desde el templado hasta semicálido, con 
precipitaciones y humedad en la estación de verano, resultando ser los 
más altos del Perú (SENAMHI, 2021). La región está particularmente 
influenciada por anomalías interanuales de El Niño, con un mayor 
transporte de humedad hacia el oriente durante El Niño, y con regímenes 
mejorados de flujo seco del occidente durante La Niña (Bergmann et al., 
2021; Perry, Seimon, 8 Kelly, 2014; Lavado €: Espinoza, 2014). 


Materiales y métodos 


Datos utilizados 


El conjunto de datos recopilado consiste en las precipitaciones diarias 
máximas anuales (PDMA), provenientes de 29 estaciones meteorológicas, 
distribuidas dentro y cerca de la región Cusco (Tabla 1). Estaciones que 
son gestionadas por el Servicio Nacional de Meteorología e Hidrología del 
Perú (SENAMHI). Algunas series de datos observados a partir de las 
estaciones presentaron falta de datos debido a periodos no operativos. 
Estos registros se sometieron a control de calidad, homogeneización y 
complementación de datos faltantes mediante el uso del paquete 
desarrollado en el software R denominado CLIMATOL (Guijarro, 2019). De 
igual forma, se recopilaron datos en formato ráster del producto IMERG 
Final Run del satélite GPM (Versión 06B), con resolución espacial de 0.12, 
para un periodo de registro de 2000 a 2019, que combina todos los 
pasivos de microondas (PMW) e infrarrojos (IR) de los satélites de la 


constelación GPM. Este producto se obtuvo de la página web de la 
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Administración Nacional de Aeronáutica y el Espacio (NASA) 
(https://disc.gsfc.nasa.gov/datasets/GPM_3IMERGDF_06). 
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Tabla 1. Valores de altitud y precipitación máxima anual en las 


estaciones meteorológicas utilizadas en el estudio. 


Núm Estación Altitud (msnm) | Precipitación (mm) 
1 Abancay 2 776 29.5 
2 Curahuasi 2741 43.6 
3 Quebrada Yanatile 1183 79.1 
4 Quillabamba 1001 59.9 
5 Ollachea 3 420 81.9 
6 Pilcopata 900 114.3 
7 Quincemil 651 179.2 
8 Caylloma 4 327 34.7 
9 Condoroma 4 287 31.4 
10 La Angostura 4 150 47.3 
11 Llally 3 985 32.9 
12 Pomacanchi 3 690 41.4 
13 Santa Rosa 3 957 45.3 
14 Santo Tomás 3 671 40.1 
15 Sicuani 3 534 39.4 
16 Tisco 4 195 46.3 
17 Yauri 3 927 35.3 
18 Anta Ancachuro 3 324 38.9 
19 Cay Cay 3117 40.6 
20 Ccateca 3 681 48.8 
21 Challabamba 2 803 48.6 
22 Cusco 3 288 33.9 
23 Granja Kayra 3 214 39.1 
24 Paruro 3 070 41.7 
25 Paucartambo 2905 54.2 
26 Pisac 2 990 35.2 
27 | Tambobamba 3 279 41.9 
28 Urcos 3 666 25.0 
29 Urubamba 2 850 41.2 
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La razón principal para la adopción del conjunto de datos IMERG es 
que el GPM dispone de series de datos con resolución de hasta 30 
minutos, particularmente adecuado para el cálculo de curvas de 


intensidad, duración y frecuencia (Kasi, Maheswaran, 8 Jarajapu, 2022). 


IMERG presenta tres productos: “Early” (cuatro horas después de 
la lectura); “Late” (14 horas después de la lectura), y “Final Run”, que es 
el resultado de un procesamiento de 3.5 meses después de la lectura y 
calibrada con información observada en tierra (Tan et al., 2019). Este 
producto se recomienda para investigaciones, dado que provee 
estimaciones más confiables (Huffman et a/l., 2019; Moazami 8 Najafi, 


2021) y fue el que se usó en esta investigación. 


Metodología 


Los datos diarios obtenidos del IMERG se sometieron a pruebas de 
identificación de datos atípicos de Grubbs (Barbato, Barini, Genta, € Levi, 
2011). En la mayoría de las estaciones se encontraron datos atípicos, en 
especial en los años 2017, 2018 y 2019, que mostraron valores altos 
(sobreestimación); en 2007 mostraron valores bajos. Esta práctica 
mejoró de forma sustancial la calidad de la información (Figura 2). Esta 


sobreestimación también la observaron Zubieta et al. (2017). 
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Figura 2. Identificación de valores atípicos en los datos provenientes 
del producto IMERG. 


En la Figura 3 se detalla la metodología seguida del enfoque mixto: 
1) con base en la información de precipitación observada se efectuó el 
análisis regional de frecuencia, donde se definen regiones homogéneas 
para finalmente construir las curvas de crecimiento por región 
homogénea; y 2) recopilación de eventos máximos de lluvia horaria del 


producto satelital IMERG, a fin de obtener los coeficientes de lluvia horaria 


“a” e invarianza de escala “n”, mediante la solución de la Ecuación (2). 
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Además, se realiza una regresión lineal entre información de precipitación 


diaria observada y datos IMERG para obtener el factor Mc. 


Metodología para la construcción de curvas IDF 
para la Región Cusco 


: Mc 
1 Control de calidad Control de calidad 
Precipitación máxima diaria Precipitación máxima diaria 
anual (pluviómetros) anual Lo 

Análisis de frecuencia de Recopilación de eventos de lluvia 
precipitaciones extremas máximos con base en la 
basado en los L-momentos precipitación horaria (IMERG) 

Regiones homogéneas Mapeo del coeficiente de lluvia 

asociadas con precipitaciones horaria (a) e invarianza de escala 
máximas en 24 horas (n) - Ecuación (13) 


| 


a Resolviendo la función ID para 
Construcción de la curva de duraciones específicas 


crecimiento Ecuación (13) 
| Curvas IDF | 
Índice de avenida 
Ecuación (12) 


Figura 3. Metodología con enfoque mixto adoptada para la construcción 


de curvas IDF, Mc=Factor de corrección. 


La estimación del factor Mc y el cálculo de los parámetros a y n se 


hicieron para cada estación meteorológica; sin embargo, el factor de 
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crecimiento se estimó por región homogénea, tal como lo plantean Qamar 
et al. (2017). 


Función intensidad-duración 


Se adoptó la propuesta basada en el concepto de índice de avenida, 
enfoque regional utilizado en diferentes investigaciones (Dalrymple, 
1960; Smith, 1989; Awadallah, 2015; Qamar et al., 2017). Con base en 


este enfoque, la precipitación viene dada por: 
h(d,T) = h(d).K(T) (1) 


Donde h(d) es el valor medio de la distribución y que depende solo 
de la duración del evento, mientras que K(T) es el factor de crecimiento 
y representa el factor de incremento, resultado de la búsqueda de eventos 
atípicos para un determinado periodo de retorno (T). Los valores de T 
empleados en este análisis varían entre 2 y 1 000 años. El valor medio de 
la lluvia h(d), también llamada función ID, se puede expresar como una 
función de la duración d (en horas) del evento extrayendo los parámetros 
a y n (WMO, 1983; OMM, 1994; Qamar et al., 2017). 


h(d) = ad” (2) 


Donde a es el coeficiente de duración de la lluvia que define la 


profundidad media de la lluvia en un intervalo de tiempo igual a una hora; 
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"n" es el exponente de la invariancia de escala que gobierna la forma de 
la curva y el grado en que depende de la duración de la precipitación. Esta 
forma de curva suele ser adecuada para las curvas 1D horarias y diarias, 


mientras que no es adecuada para los extremos subhorarios. 


De acuerdo con la metodología propuesta por Qamar et al. (2017), 
para obtener la función intensidad-duración (ID) se debe aprovechar la 
información satelital, ya sea de resolución diaria, así como horaria, los 
cuales fueron posibles en esta investigación, dado que se utilizó la 
información ofrecida por el producto satelital IMERG VO6 (Huffman et al., 
2019; Moazami € Najafi, 2021). 


Para llevar a cabo el proceso de corrección de los datos satelitales 
del producto IMERG con respecto a los datos observados, se acumularon 
los datos del pluviómetro y satelitales disponibles a paso diario en 
intervalos de 1, 2, 3 y 7 días, y se ajustaron mediante una regresión 
lineal, para obtener el valor del factor de corrección o factor de sesgo Mc 
que se utiliza en la corrección de datos satelitales (Hunink, Immerzeel, € 
Droogers, 2014). 


Luego se estimaron los factores de corrección para las 29 
estaciones, obteniendo valores que varían desde 0.61 para la estación de 
Urubamba hastas 1.1 para la estación de Tambobamba. Dichos valores 
de Mc se usan como un coeficiente multiplicativo para obtener los datos 
finalmente corregidos de precipitaciones máximas diarias anuales 
provenientes del satélite en cualquier punto requerido de la región Cusco. 
Los valores de Mc muestran que los datos del IMERG presentan valores 


superiores a los observados. 
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Por último, con el uso del paquete R y con base en la información 
descargada del IMERG de precipitaciones máximas con distinta duración, 
se pudieron seleccionar para cada estación los eventos de lluvia para 
efectuar las regresiones que conduzcan a la obtención de los parámetros 
a y n, que devienen de la linealización de la Ecuación (2), resultando el 
espacio log h vs log d para duraciones de 3, 6, 9, 12 y 24 horas, y así 


contar con las curvas de intensidad y duración (ID). 


Análisis regional de frecuencia y regiones homogéneas 


Para identificar áreas con similar comportamiento hidrológico (regiones 
homogéneas), se siguió el procedimiento recomendado por Hosking y 
Wallis (1997), que consiste en un agrupamiento inicial y refinamiento 
manual de regiones. Para dicho propósito se emplearon predictores físicos 
estandarizados (altitud, latitud, longitud y precipitación media anual); 
luego se adaptó el método de agrupamiento jerárquico de Ward en 
consideración con las estadísticas de los L-momentos (Charrad, Ghazzali, 
Boiteau, 8: Niknafs, 2014). Finalmente, basado en el último agrupamiento, 
se acondicionan las regiones considerando los estadísticos de discordancia 
(Di) y heterogeneidad (H). 


Los L-momentos introducidos por Hosking (1990) se utilizan 
ampliamente para ajustar funciones de distribución al análisis regional 
(Stedinger, Vogel, Foufoula-Georgiou, 1993; Hosking 8 Wallis, 1997; 
Anctil, Larouche, Viau, € Parent, 2002; Yurekli, 2005). La linealidad de 


los L-momentos, en comparación con los momentos convencionales, tiene 
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la ventaja de ser insensible a los valores extremos (Stedinger et a/., 1993; 
Sankarasubramanian € Srinivasan, 1999). 


Los L-momentos son combinaciones lineales de momentos 
ponderados (MP), según la definición de Hosking (1990), y Landwehr, 
Greenwood, Matalas y Wallis (1979), la cual resume la teoría de los 


momentos ponderados de la siguiente manera: 
BP, = EfX[E.(0 1" (3) 


Donde 6B, es el is, orden de los MP; EfX), el valor esperado; Fx(X), la 
función de distribución acumulativa, y r es el valor r-ésimo. Hosking y 
Wallis (1997) definieron un estimador insesgado de una muestra MP 
designada (b;) y obtuvieron un estimador insesgado de la muestra de los 
primeros cuatro L-momentos por el estimador de la muestra de MP. La 
muestra no sesgada estimada de MP para cualquier distribución se puede 


calcular a partir de Landwehr et a/. (1979): 


Bn == x; (4) 


r 


n- l 


Donde (x;) es un conjunto ordenado de observaciones como X: <€ X2 
< X3< ... Xn, donde n es el número de datos de muestra y j es el valor j- 
ésimo. Para cualquier distribución, los primeros cuatro L-momentos se 


calculan con facilidad a partir de los MP: 
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A1= bi 
A2= 2b, — b; 
23= 6b3 — 6b, + b, 
»24= 20b4 — 30b3 + 12b, — b, (5) 


Donde xa, es la media de la distribución; »,, la medida de escala; 
13, la medida de desviación, y »4 es la medida de curtosis. Las relaciones 
de los L-momentos (L-coeficiente de variación, L-asimetría y L-curtosis) 


se definen como sigue: 


E (6) 


Para validar la homogeneidad de una región en términos de 
relaciones de L- momentos se utilizó la prueba de discordancia propuesta 
por Hosking y Wallis (1993). 


p E ay T . 
Sea u = |t0,1P,e7?] el vector que contiene los vectores t, ts, ta en 
el sitio /, donde el superíndice T denota transposición de un vector o 


matriz. 


Sea u=N"*Y*,u, el promedio regional no ponderado de los L- 
momentos para cada grupo, donde N es el tamaño de la muestra en cada 


grupo. La medida de discordancia para el sitio ¡ se define entonces como: 
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D; =5(u; - "Su; —1) (7) 


En la ecuación anterior, S”* es la inversa de la matriz S, donde S = 
WD, -0D(, — 34). Hosking y Wallis (1997) propusieron el criterio 
de D; > 3 (Ecuación (7)) para excluir una estación de la región 
homogénea. La relación entre sitios se utiliza para identificar sitios 
coherentes con una distribución de frecuencia similar. La prueba de 
heterogeneidad H, compara las muestras de las relaciones de los L- 
momentos con los parámetros de distribución kappa; mide la 
heterogeneidad entre sitios en la misma región. Hosking y Wallis (1997) 


propusieron la siguiente estadística: 


1 =*2 (8) 


Op 


Donde V y gd, son la media y desviación estándar de los valores 


simulados de V,. Los V, se calculan mediante la siguiente ecuación: 


v, = (E, m [e 0/12," (9) 


Las simulaciones se llevan a cabo utilizando una distribución flexible 


con el promedio regional de las relaciones de los L-momentos 1, T, T3 y 
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Ta. Siguiendo a Hosking y Wallis (1997), se usó la distribución kappa de 


cuatro parámetros con la función de cuantiles: 
x(F) = E +1 —[(1 — F*)/A1)/k (10) 


Medidas de heterogeneidad, como la estadística H1, han sido 
desarrolladas por Hosking y Wallis (1997), como un indicador del grado 
de heterogeneidad en los L-momentos para un grupo de sitios/estaciones. 
La estadística H1 mide la variabilidad relativa del coeficiente de variación 
(L-CV), observado a partir de la muestra, y se utiliza para analizar la 
heterogeneidad en cada una de las regiones homogéneas propuestas 
(Wallis, Schaefer, Barker, € Taylor, 2007). 


Existen dos medidas adicionales de heterogeneidad, denominadas 
H2 y H3. Sin embargo, no se recomienda su uso debido a que se ha 
comprobado que presentan una notable falta de potencia estadística 
(Viglione et al., 2007). El fundamento teórico y procedimiento de cálculo 
de la medida de heterogeneidad H1 está ampliamente descrito por 
Hosking y Wallis (1986). 


La región se considera “aceptablemente homogénea” si H < 1; 
“posiblemente heterogénea” si 1 < H < 2, y "definitivamente 


heterogénea” si H > 2: 
Y H1: medida de homogeneidad en términos de L-CV, 
Y H2: medida de homogeneidad en términos de L-CS. 


Y H3: medida de homogeneidad en términos de L-CK. 
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El estadístico |ZPIST| juzga qué tan bien la L-asimetría simulada y la 
L-curtosis de una distribución ajustada coinciden con el promedio regional 
de L-asimetría y L-curtosis que fueron obtenidos de valores de los datos 
observados (Hosking, 1990). El criterio de bondad de ajuste para cada 
distribución se define mediante el estadístico |zP”!"|, en función de las 


distintas distribuciones candidatas: 
ZST = (1% — 14 + Ba)/0, (11) 


Con tf = valor promedio de ta de los datos de la región Ba; 04 = 
sesgo y desviación estándar de ra, respectivamente, y se definen de la 


siguiente manera: 


By=Nh = Er” 4) (12) 
- -1 Nsim( (mM) R ; 2 iS 
04 = WNsim E 1) ala 7 1) S Nsim Ba (13) 


Donde Nim = número fijo de simulaciones de datos regionales 
generados utilizando una distribución Kappa (Hosking, 1990); m = mth, 
región simulada obtenida utilizando una distribución Kappa. El ajuste se 
considera adecuado si |ZP!5T| está lo suficiente cerca de cero. Si |ZP!5T| < 
1.64, el ajuste es razonable. Si es aceptable más de una distribución 


candidata, la que tenga el |ZP'5'| más bajo se considera la más adecuada. 
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Además, la elección de la función de distribución y la bondad de ajuste de 
las funciones de distribución se verifican utilizando el diagrama de L- 
momentos (Hosking 8 Wallis, 1997), para cada grupo derivado de la 
regionalización de la precipitación diaria máxima anual de las estaciones 
en toda el área de estudio. La idea detrás del uso de diagramas de L- 
momentos se basa en el aprovechamiento de combinaciones únicas de 
coeficientes de asimetría y curtosis, a fin de identificar de forma gráfica 
la función más cercana a la muestra de estudio (Kumar, Chatterjee, 
Kumar, Lohani, € Singh, 2003; Lin € Chen, 2006). 


Después de haber identificado la distribución de frecuencia 
adecuada para cada región homogénea se estiman los cuantiles para 
diferentes periodos de retorno utilizando el método del índice de avenida. 
El método asume que los datos de precipitación máxima diaria en 
diferentes sitios en una región homogénea presentan la misma 
distribución, excepto por un parámetro de escala específico del sitio o un 
factor índice (Dalrymple, 1960). El factor de escala se denomina como un 
índice de avenida y por lo general se toma como la media de la 
precipitación diaria máxima anual. Este hecho hace que la información 
resultante con base en esta metodología sea útil para zonas donde no 


existen estaciones meteorológicas. 


Las estimaciones de cuantiles Ó(F) con probabilidad de no 
excedencia F en un sitio en una región con N sitios se calculan mediante 
Q((P) =1K(M(P), donde K(T) es una función cuantílica adimensional 
común (curva de crecimiento regional), y "!,” es el valor índice de avenida, 
que representa el cuantil T-año de la distribución regional normalizada 
(Ngongondo et a/., 2011). 
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La disponibilidad de datos de precipitación del producto satelital 
IMERG (VO6B) permitió obtener una gran variedad de registros de 
tormentas que varían en intensidad y duración. Estos registros permiten 
analizar eventos de lluvia extrema y finalmente estimar los parámetros 
regionales a (coeficiente de duración de lluvia) y n (coeficiente de 
invarianza escalar); parámetros útiles para optimizar tanto la profundidad 


de la lluvia como la forma de la curva. 


Empleando los parámetros regionales para la distribución 
identificada, la curva de crecimiento regional se calcula y se multiplica por 
la precipitación máxima anual promedio específica de la estación para 
obtener los cuantiles de precipitación deseados para la estación relevante. 
Para cada región homogénea se simula una región que tiene el mismo 
número de estaciones, registro de datos en cada estación, 
heterogeneidad, promedio regional y relaciones de L-momentos promedio 
regionales como los datos observados. Asimismo, para cada simulación 
se calcularon los errores en la curva de crecimiento simulada y los 
cuantiles, y luego se estimaron el sesgo, error cuadrático medio (RMSE) 


y límites de error del 90 %. 


Previo a analizar la relación entre datos de lluvia observada y 
satelital se empleó el factor de corrección Mc para mejorar la calidad de 
los datos estimados provistos por el producto IMERG. Finalmente, se 
realizó una regresión lineal entre los máximos anuales basados en datos 


observados y satelitales para la duración diaria (Qamar et al., 2017). 
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Resultados 


Regiones homogéneas 


Luego del análisis regional de frecuencia basado en los L-momentos se 
definieron cuatro regiones bajo los siguientes criterios: una región se 
considera “aceptablemente homogénea” si Hi < 1; “posiblemente 
heterogénea” si 1 < H:¿< 2, y “definitivamente heterogénea” si Hi > 2 
(Tabla 2 y Figura 4). Quiere decir que la hipótesis de agrupamiento 
efectuada es aceptablemente homogénea para las cuatro regiones dado 
que en la columna H; de la Tabla 2, los valores para R1 a R4 son menores 
al. 


Tabla 2. Valores de H1, H2 y H3 para las regiones homogéneas. 


Regiones Hi H2 H3 
R1 0.1508 -1.2183 0.1543 
R2 0.2250 1.5889 1.5057 
R3 -0.0324 -0.3788 0.1326 
R4 0.4967 1.1569 2.7210 


Nota: Hi H2 y H3 se refieren a la medida de homogeneidad en términos de L- 


coeficiente de variación, L-asimetría y L-curtosis, respectivamente. 
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Altitud Precipitación 


Nro Estación Abreviatura 
(msnm) (mm) 
1 A Abancay ABAN 2776 29.5 
2 A  Curahuasi CURA 2741 43.6 
3 4 Quebrada Yanatile QUEB 1183 79.1 
4 A Quillabamba QUIL 1001 59.5 
5 O  Ollachea OLLA 3420 81.9 
6 e  Pilcopata PILC 900 114.3 
ax 7 e  Quincemil QUIN 651 179.2 
” 8 m  Caylloma CAYL 4327 34.7 
MADRE DE DIOS 9 Mm  Condoroma COND 4287 31.4 
10 Mm La Angostura ANGO 4150 47.3 
11 m  Llally LLAL 3985 32.9 
12 m  Pomacanchi POMA 3690 41.4 
13 M. Santa Rosa li SANT 3957 45.3 
14 M- Santo Tomas SANT 3671 40.1 
15 Mm Sicuani SICcUu 3534 39.4 
e 16 m  Tisco TISC 4195 46.3 
* 17 E - Yauri YAUR 3927 35.3 
18 * Anta Ancachuro ANTA 3324 38.9 
19 * Cay Cay CAY 3117 40.6 
20 *+  Ccatcca CCAT 3681 48.8 
21 x  Challabamba CHAL 2803 48.6 
22 * Cusco cusc 3288 33.9 
23 * Granja Kcayra GRAJ 3214 391 
APURIMAC S 24 x Paruro PARU 3070 41.7 
Z 25 x  Paucartambo PAUC 2905 54.2 
26 *  Pisac PISA 2990 35.2 
27 k*  Tambobamba TAMB 3279 41.9 
28 *  Urcos URCO 3666 25 
29 * Urubamba URUB 2850 41.2 
Ay LEYENDA 

AYACUCHO e la CA Región Cusco 

SE Regiones homogeneas 
kilometros 
AB 0 E  x 
1 2 3 4 


-73* -72 


Figura 4. Distribución espacial de las estaciones meteorológicas por 


región homogénea. 


Con base en el estadístico |ZPIST| se describe que la función de 


distribución de probabilidad GLO se ajusta a las regiones 1 y 2, mientras 


la distribución GNO se ajusta a las regiones 3 y 4 (Tabla 3). 
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Tabla 3. Valores de [ZP*!9"] para cada distribución evaluada por región 
homogénea. 
[ZP157] 
p Valor 
Logístico Normal Pearson Pareto . 
0 ] extremo p . ] Mejor 
Región | generalizado ] generalizado tipo III | generalizado 
generalizado FDP 
(GLO) (GNO) (PE3) (GPA) 
(GEV) 
Ri -1.5792 -2.7285 -2.8571 -3.2205 -5.2662 GLO 
R2 -0.7625 -1.4159 1,9995 -1.9636 -2.0551 GLO 
R3 1.7853 -0.5513 -0.5082 -0.8762 -5.4179 GNO 
R4 1.8727 -1.3917 -0.8087 -0.8832 -7.6164 GNO 


Según lo planteado por Hosking y Wallis (1997), se considera que 


el ajuste de una distribución es aceptable si el valor del estadístico |ZP15"] 


es cercano a cero. La estimación de los cuantiles se realiza una vez 


escogida una función de distribución, aplicando la prueba de bondad de 


mejor ajuste. 


En la Figura 5 se aprecia la ubicación de los L-momentos promedio 


regionales con ratios teóricos de L-asimetria/L-curtosis para las diferentes 


FDP. De esta manera se corrobora que las FDP ubicadas en las regiones 


1 y 2 están más cerca de la curva GLO y GEV, por lo que la distribución 


GLO resulta con mejor estadístico [ZP!57|, En el caso de las regiones 3 y 4, 


las curvas GNO y PE3 presentan mejores ajustes; sin embargo, GNO fue 


superior. 
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0,4 


0.2 e 
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Figura 5. Diagrama de los L-momentos de las distribuciones en las 
regiones identificadas: a) región 1, distribución elegida GLO; b) región 
2, distribución elegida GLO; c) región 3, distribución elegida GNO, y d) 

región 4, distribución elegida GNO. 


Se determinaron los factores de crecimiento (FC) para los periodos 
de retorno (Tr) de 2, 5, 10, 20, 50, 100, 200, 500 y 1 000 años (Tabla 4 
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y Figura 6). Los resultados sugieren una mayor incertidumbre en la región 
2; de hecho, para un periodo de retorno de 1 000 años se tienen valores 
de FC más altos que en otras regiones. Esto se debe tal vez a la limitada 
serie de datos observados. En contraste, la región 4 presenta los menores 
valores de FC. Asimismo, se puede observar que la región 1 tiene un 
comportamiento similar a la región 2; la región 3 presenta condiciones 
intermedias. Esta similitud entre las regiones 1 y 2 pueden deberse al 
régimen climático estacional, ya que se ubican sobre todo en zonas 
amazónicas; mientras que las regiones 3 y 4 están ubicadas en zonas 


andinas. 


Tabla 4. Valores del factor de crecimiento (Tr: periodo de retorno). 


Tr (años) Región 1 Región 2 Región 3 Región 4 
2 0.96 0.95 0.97 0.99 
5 1.18 1.17 1.19 1.20 
10 1.34 1.33 1.32 1.33 
20 1.51 1.50 1.44 1.43 
50 1.75 1.77 1.59 1.55 
100 1.95 2.00 1.70 1.64 
200 2.18 2.27 1.81 1.71 
500 2.53 2.69 1.94 1.81 
1000 2.82 3.07 2.05 1.88 
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Curva de crecimiento 


10 20 50 100 200 500 1000 2 5 10 20 50 100 200 
Periodo de retorno (Tr) Periodo de retorno (Tr) 


Curva de crecimiento 


10 20 50 100 200 500 1000 2 5 10 20 50 100 200 


Periodo de retorno (Tr) Periodo de retorno (Tr) 


Figura 6. Curva de crecimiento adimensional por región homogénea. 


Relación entre datos de lluvia tanto satelital como 
observados y factor Mc 


Para analizar el vínculo entre la lluvia observada y aquella obtenida por el 
satélite IMERG, se seleccionaron dos estaciones (Urubamba y 
Tambobamba). Los valores encontrados para Mc variaron desde 0.61 para 


la estación de Urubamba hasta un valor de Mc = 1.1 para la estación de 
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Tambobamba (Figura 7). Tales resultados muestran la importancia de la 
aplicación del factor de corrección dado, que valores como 0.61 sugieren 


la necesidad de corregir información satelital para su aplicación. 


¿AA A A OA OA EXA Ez? 2X¿»AE5EE EA ——— 


A 


os | o_s 
o ] 
ESTACION TAMBOBAMBA e % " [ ESTACION URUBAMBA Pr ] 
96 | 150! 
Ml 1 dia M 1dia 
A 2dias | 1% Yo A 2dias | 
A “Y 3 dias ; | Y y $ / 3dias 
dk e a O 7dias | | O O 7 dias | 
rá ya | 50)! Ma A 
UE 4 ñ R3==0.98 | | ue = R? == 0.99 
E me Mc == 1.1 [ A M6 == 0.61 
ET INEA EE SI A ON E O 0 aci ñ A TU A A A A E A 
50 100 150 200 250 300 0 50 100 150 200 250 300 
Precipitación anual máxima IMERG a diferentes intervalos (mm) Precipitación anual máxima IMERG a diferentes intervalos (mm) 


Figura 7. Estimación del factor de corrección Mc con base en los datos 
diarios observados, y datos IMERG para las estaciones de Tambobamba 


y Urubamba. 


A diferentes intervalos (1, 2, 3, 7 días) se observa una relación 
significativa entre la precipitación anual máxima obtenida a partir de 
datos satelitales IMERG y datos observados, para las estaciones 
Tambobamba (R? « 0.98, p < 0.01) y Urubamba (R? « 0.99, p < 0.01), 


estaciones ubicadas en el centro del área de estudio (Figura 1). 
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Distribución espacial de eventos extremos: función 
intensidad-duración 


Inicialmente se estimó el coeficiente de duración de lluvia a, el cual 
presenta valores menores a 40 (considerados valores bajos) en la zona 
andina de la región del Cusco; mientras que en la selva de Quincemil, 
Pilcopata y Tambobamba presenta valores mayores a 60 (considerados 
valores altos). Por otro lado, el coeficiente de invarianza escalar n 
presenta valores menores a 0.15 (considerados valores bajos) en la parte 
oeste de la región, y valores mayores a 0.15 en la parte norte y sur este 


de la región (Figura 8). 
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Leyenda Leyenda 
3 Región Cusco C3 Región Cusco 
Regiones Regiones 
homogéneas homogéneas 
A  R1:Región 1 A  R1:Región 1 
e  R2:Región 2 e  R2:Región 2 
n R3 : Región 3 n R3 : Región 3 
* R4:Región 4 *  R4: Región 4 
74" 73 72" 71? 74? TEN 72* ri 
Longitud Longitud 


Figura 8. Mapeo del coeficiente de duración de lluvia horaria "a" y del 


exponente de la invariancia de escala "n" en la región Cusco. 


Los valores de a encontrados para la zona de sierra de la región del 
Cusco son bajos, menores a 40 (la zona andina del Cusco registra una 
precipitación máxima anual media de 45 mm/día), mientras que en la 


zona amazónica del Cusco presenta valores altos mayores a 60 (la selva 
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del Cusco registra una precipitación máxima anual media de 100 mm/día) 
(Figura 8). 


Bajo la hipótesis de que el valor del factor Mc se mantiene constante 
para diferentes duraciones de tiempo, el valor medio de la lluvia (Ecuación 
(2)) h(d) puede ser multiplicado por el factor Mc para obtener la 
aproximación de la lluvia en el punto de interés requerido. Así, 
posteriormente fue posible obtener la distribución espacial de h(d) para la 


región Cusco (Figura 9). 
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Figura 9. Mapeo de h(d) en la región Cusco para duraciones de 3, 6, 12 


y 24 horas. 
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Con base en estos valores obtenidos, se generan mapas ráster 
mediante interpolación para periodos de retorno de 20 y 100 años, y para 
duraciones de 3, 6, 12 y 24 horas (Figura 10 y Figura 11). El enfoque 
mixto (que emplea datos IMERG en regiones donde no se cuenta con 
datos observados) propuesto en esta investigación muestra mapas cuyos 
resultados se aproximan más a los datos observados, validando de esta 
forma la metodología propuesta. En los mapas se aprecian valores altos 
en las zonas de Quincemil y Chontachaca; valores bajos en las localidades 
de Espinar y Cusco; y valores medios y transiciones en la ceja de selva, 


selva norte y sectores próximos a Challabamba y Paucartambo. 
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Figura 10. Mapas de lluvias extremas para duraciones de 3, 6, 12, y 24 


horas y periodo de retorno de 20 años. 
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Figura 11. Mapas de lluvias extremas para duraciones de 3, 6, 12, y 24 


horas y periodo de retorno de 100 años. 
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Incertidumbre de los resultados obtenidos 


Se compararon valores a partir del uso del enfoque mixto y del producto 
IMERG con valores obtenidos desde pluviómetros. En contexto, se 
encontró que el enfoque mixto presenta valores del RMSE más cercanos 
a cero y valores de RMSE menores a los obtenidos a partir del producto 
IMERG. Asimismo, las estaciones que presentan valores más bajos de 
RMSE son Granja Kayra, Paruro y Tambobamba, en la región central oeste 
del Cusco (región 1). En contraste, la región 4, donde se encuentra la 


estación de Quincemil, presenta valores más altos (Tabla 5). 


Tabla 5. Valores del error cuadrático medio (RMSE) para diferentes 


periodos de retorno en las estaciones representativas de cada región. 


Región Estación Tr Año IMERG Enfoque mixto 
2 3.19 1,04 
5 4,09 1.27 
10 4,64 1.40 
20 5.19 1.51 
Granja Kayra 50 5.97 1.64 
100 6.62 1.72 
200 7.34 1.81 
500 8.37 1.91 
Región 1 1 000 9.23 1,99 
2 3.58 0.64 
5 4.67 0.78 
10 5.36 0.85 
20 6.05 0.92 
Urubamba 50 7.02 1.00 
100 7.82 1.05 
200 8.67 1.10 
500 9.89 £.12 
1000 10.88 1.21 
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Región Estación Tr Año IMERG Enfoque mixto 
2 3.02 1.22 
5 3.85 1,49 
10 4.42 1.64 
20 5.02 1.77 
Paruro 50 5.91 1.92 
100 6.67 2.03 
200 7.49 2.13 
500 8.69 2.25 
1 000 9.68 2.34 
2 3.22 0.58 
5 4,13 0.70 
10 4.74 0.78 
20 5.37 0.84 
Tambobamba an jóS 9197 
100 7.04 0.96 
200 7.85 1.00 
500 9.02 1.06 
1 000 9.96 1.10 
2 4.48 0.66 
5 6.11 0.81 
10 7,44 0.92 
20 9.06 1.03 
Quillabamba 50 11.87 1.20 
100 14.67 1,34 
200 18.18 1,49 
500 24.12 1.73 
1 000 29.77 1.93 
REGION 2 2 6.67 0.71 
5 8.41 0.88 
10 9.89 1.00 
20 11.76 1.12 
Quebrada 50 15.13 1.30 
Yanatile 100 18.58 1.45 
200 22.95 1.62 
500 30.43 1.88 
1 000 37.58 2.10 
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Región Estación Tr Año IMERG Enfoque mixto 
2 4.30 0.24 
5 5.66 0.29 
10 6.49 0.33 
20 7.29 0.36 
Santo Tomas cd 209 gan 
100 9.27 0.42 
200 10.20 0.45 
500 11.51 0.48 
1 000 12.57 0.51 
2 4.30 0.24 
5 5.66 0.29 
10 6.49 0.33 
20 7.29 0.36 
Caylloma 50 8.39 0.40 
100 9.27 0.42 
200 10.20 0.45 
500 11.51 0.48 
1000 12.57 0.51 
Regidh:a 2 3.65 2.82 
5 4.82 3,45 
10 5.62 3.84 
20 6.45 4,19 
La Angostura 50 7.66 4.63 
100 8.65 4.94 
200 9.71 5.25 
500 11.21 5.66 
1 000 12.42 5.96 
2 3.92 0.98 
5 5.11 1,19 
10 5.93 1.33 
20 6.78 1.45 
a 50 8.01 1.60 
100 9.03 1.71 
200 10.12 1.82 
500 11.67 1.96 
1 000 12.91 2.06 
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Región Estación Tr Año IMERG Enfoque mixto 
2 6.16 0.29 
5 6.60 0.36 
10 7.12 0.41 
20 7.90 0.46 
Ollachea 50 9.45 0.54 
100 11.09 0.61 
200 13.16 0.69 
500 16.56 0.82 
Región 4 1 000 19.64 0.94 
2 17.79 1.34 
5 20.31 1.64 
10 22.89 1.87 
20 26.06 2.11 
Quincemil 50 31.33 2.49 
100 36.25 2.81 
200 42.03 3.19 
500 51.07 3.78 
1 000 59.04 4.31 


Por otro lado, la metodología del estudio se evaluó mediante el 
método de validación cruzada entre los valores de precipitación diaria 
obtenidos del producto IMERG, enfoque mixto y precipitación observada 
(Figura 12). Esto se realizó para periodos de retorno de 2, 5, 10, 20, 50, 
100, 200, 500 y 1 000 años; se puede notar en la mayoría de los casos 
que existe sobreestimación en los valores extremos de lluvia al emplear 
datos provenientes del producto IMERG. Esta sobreestimación de valores 
extremos es más notoria para las estaciones de Urubamba, Quillabamba, 
Quincemil y Quebrada Yanatile (a excepción de la estación Tambobamba, 
que presenta data similar a la observada). Mientras que los valores que 
se obtienen mediante el enfoque mixto resultan muy parecidos a la data 


observada. 
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Figura 12. Comparación de valores de precipitación obtenidos mediante 


datos del IMERG, enfoque mixto y datos observados para diferentes 


periodos de retorno y estaciones representativas. 


Discusión 


La densidad promedio de estaciones consideradas en esta investigación 


hace que el área de influencia por estación meteorológica sea de 2 482 


km, 


Organización Meteorológica Mundial 
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considerada muy superior a la densidad recomendada por la 
(OMM), 


la cual indica que las 
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distancias entre estaciones debieran ser de 25 a 30 km para zonas planas 
y de 12.5 a 15 km en zonas de montaña (densidad promedio de 225 km? 
por estación) (WMO, 1970; OMM, 2010). Ante ello, resulta relevante 
incrementar el número de estaciones meteorológicas en el departamento 
del Cusco y del Perú (Rojas et al., 2021). No obstante, en las últimas 
décadas surge la importancia y/o necesidad de evaluar la aplicabilidad de 
datos satelitales en la estimación de eventos hidroclimáticos extremos en 
zonas andinas y amazónicas del Perú con fines de prevención (Zubieta et 
al., 2019). Tal es el caso del producto en IMERG (Final Run, GPM Versión 
06B) del satélite GPM, el cual es recomendado para investigaciones, dado 
que provee datos más adecuadas para la estimación de variables 
hidrológicas (Zubieta et al., 2017; Huffman et al., 2019; Moazami 8 
Najafi, 2021). 


Para estimar la distribución de probabilidad de lluvias extremas ante 
diferente duración (días) y diferente tiempo de retorno. Nuestro estudio 
sugiere que para alcanzar un mayor rendimiento los valores de 
precipitación satelital deben ser sometidos previamente a un factor 
multiplicativo (Mc). La estimación de este factor Mc y el cálculo de los 
parámetros necesarios a y n se hizo por cada estación meteorológica; sin 
embargo, el factor de crecimiento se valoró por región homogénea. Este 
estudio pone a disposición mapas de lluvias máximas extremas (duración 
de 3, 6, 12 y 24 horas), correspondiente a periodos de retorno de 20 y 
100 años. Dichos mapas pueden ser útiles para el diseño de 
infraestructura hidráulica para sistemas de drenaje pluvial o control de 


inundaciones (Fernández 8: Lavado, 2017). 
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El empleo de productos IMERG para análisis de frecuencias y 

representación de caudales extremos puede conducir a rendimientos 
satisfactorios (Jiang et al., 2022; Chen, Zhu, Wang, 8l Liao, 2022). De 
hecho, en nuestro estudio, en comparación con datos observados, el 
análisis de RMSE muestra que los valores de precipitación obtenidos a 
partir del enfoque mixto presentan valores más adecuados para 
representar lluvia extrema, que aquellos obtenidos mediante el empleo 
de solo datos satelitales IMERG. 


Las regiones homogéneas 1 y 2, ubicadas en la zona central y sur 
del Cusco, presentan los valores más bajos de RMSE. Ello sugiere que el 
uso del enfoque mixto es más adecuado para la representación de eventos 
de lluvia extrema bajo diferentes tiempos de duración y diferente tiempo 
de retorno en zonas andinas. Esto se asocia posiblemente con eventos de 
convección proveniente de flujos de humedad del este, muy frecuente 
durante la temporada de lluvias en zonas de alta montaña (Flores-Rojas 
et al., 2021). Esto es consistente con la fuerte estacionalidad de la lluvia 
en zonas andinas, la cual se rige en la región de los Andes peruanos 
(Espinoza et al., 2009). 


Del mismo modo, se aprecia que las estaciones de la región 1 y 2 
presentan mejores valores de RMSE que las regiones 3 y 4. En 
comparación a las regiones 1 y 2, la reducida representación de eventos 
de lluvia extrema (mayores valores de RMSE) en las regiones 3 y 4 sugiere 
que el uso del enfoque mixto sobre las regiones amazónicas puede ser 
limitado. Tal nivel de incertidumbre se puede deber a la limitada cantidad 
de estaciones meteorológicas analizadas sobre regiones amazónicas en el 


este del Cusco. De hecho, en contexto, la disponibilidad de datos 
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climáticos para investigaciones en el Perú se caracteriza por contar con 
mayor número de estaciones principalmente en zonas andinas que en 


zonas amazónicas (Zubieta et a/., 2017; Lavado et a/., 2012). 


La comparación del número de eventos extremos estimados a partir 
de datos IMERG sugiere una sobreestimación con respecto a datos 
observados en la cantidad de eventos extremos identificados, en 
particular en zonas amazónicas del norte y este del Cusco. Ello es 
consistente con la distribución espacial de lluvias extremas para 
diferentes duraciones (3, 6, 12 y 24) y diferentes periodos de retorno (20 
y 100 años), que muestra que las más intensas lluvias, como era 
esperado, es sobre regiones amazónicas. No obstante, el uso del enfoque 
mixto dispone una mejor aproximación de estos eventos extremos de 
lluvia, ya sea en zonas andinas o amazónicas. Nuestros hallazgos sugieren 
que el enfoque mixto, como método alternativo, puede utilizarse para 
generar información de IDF en zonas no monitoreadas o carente de datos 
(Fernández 8i Lavado, 2017). 


Conclusiones 


Esta investigación documenta la aplicación del análisis regional de 
frecuencias de precipitaciones máximas en 24 horas; se proponen 4 
regiones homogéneas para la región Cusco ubicada en los Andes del sur 
del Perú. Basados en series de tiempo de lluvia (observada y satelital 
IMERG-GPM) y predictores físicos estandarizados, como altitud, latitud y 
longitud, se empleó la técnica L-momentos para finalmente estimar 


funciones de distribución de probabilidad con el mejor ajuste para cada 
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región homogénea, que condujeron a la estimación de distribución 
espacial de eventos de lluvia extrema para diferentes duraciones y 
diferente tiempo de retorno. A partir de la precipitación de satélite, las 
regiones homogéneas posibilitan estimar valores de intensidad duración 


y frecuencia. 


La resolución temporal de datos satelitales IMERG (30 minutos) 
permite identificar eventos de tormentas, que junto con los datos 
observados contribuyó a proponer un método de enfoque mixto 
(observado y satelital), para la estimación de lluvias extremas en zonas 


andino y amazónicas de la región Cusco en el sur del Perú. 


En comparación con datos observados, nuestros hallazgos indican 
que el enfoque mixto provee resultados más adecuados que los datos 
satelitales IMERG. Esto sugiere que la metodología adoptada es una 
propuesta alternativa para estimar valores de máximas precipitaciones en 
lugares no monitoreados o carentes de información de la región Cusco. 
Los resultados más precisos pueden ser alcanzados en zonas andinas que 


en amazónicas. 


El análisis de incertidumbre desarrollado sugiere que el uso de los 
resultados de este estudio con fines de diseño de ingeniería debe ser 
tomado considerando las bandas de confianza de las curvas de 
crecimiento regional, además de prever que el periodo de las series de 
tiempo utilizadas en algunas regiones es corto, y que el número de 
estaciones para el análisis regional de frecuencia fue muy limitado para 


las regiones amazónicas del Cusco. 
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